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James  K.  Freericks 
November,  29,  2005 


Summary: 

This  is  the  final  technical  report  for  James  Freericks'  ONR  superconducting  electronics 
grant  for  the  period  from  1999—2005.  This  work  involved  the  development  of  computer 
modeling  codes  to  describe  Josephson  junctions  with  barrier  materials  that  are  tuned 
close  to  the  metal-insulator  transition.  The  figure  of  merit  can  be  determined  as  a 
function  of  different  parameters,  including  the  proximity  to  the  metal-insulator  transition, 
the  thickness  of  the  barrier,  and  the  temperature.  We  discovered  a  simple 
phenomenological  parameter,  called  the  generalized  Thouless  energy,  that  is  useful  for 
characterizing  and  summarizing  the  Josephson  junction  properties;  much  behavior  of 
these  systems  becomes  universal  when  represented  in  terms  of  this  Thouless  energy.  We 
also  developed  a  nonequilibrium  many-body  formalism  that  can  determine  the  response 
of  these  systems  to  large  electric  fields.  Much  of  the  computational  work  was  carried  out 
on  DOD  parallel  machines  from  the  HPCMO  facilities. 

Summary  of  technical  projects  completed 

A.  Modeling  of  Josephson  junctions. 

The  major  technical  project  completed  in  this  grant  was  the  development  of  a  computer 
modeling  program  that  can  describe  the  behavior  of  Josephson  junctions  with  the  barriers 
tuned  to  lie  close  to  a  metal-insulator  transition.  This  regime  is  of  interest  for  self- 
shunted  high  speed  JJs  that  would  be  used  for  ultrafast  digital  electronics.  The  main 
experimental  system  that  was  being  worked  on  was  materials  with  TaxN  barriers,  where 
the  material  undergoes  and  MIT  when  x=0.6.  Our  approach  was  based  on  inhomogeneous 
dynamical  mean  field  theory,  which  is  a  nonperturbative  technique  that  can  determine  the 
many-body  Green’s  functions  and  transport  of  multilayered  nanostructures.  We 
developed  a  series  of  calculations  that  can  describe  the  superconducting  state  and  the 
situation  where  the  superconductor  carries  current  due  to  a  phase  gradient  placed  over  the 
superconducting  leads.  This  allows  for  a  calculation  of  the  current  density  of  the  JJ. 

Next,  we  developed  a  real-frequency,  normal-state  code,  that  can  calculate  the 
conductivity  matrix  via  a  Kubo  formula,  from  which  the  resistance-area  product  can  be 
determined.  The  product  of  these  two  quantities  yields  the  JJ  figure  of  merit.  We 
examined  a  number  of  different  cases  through  our  studies  including:  (i)  tuning  the  barrier 
through  a  MIT  by  varying  the  Coulomb  interaction;  (ii)  examining  the  effects  of 
electronic  charge  reconstruction  that  can  drive  the  junction  into  a  superconductor- 
insulator-normal  metal-insulator-superconductor  (SINIS)  structure;  (iii)  temperature 
dependence  of  the  figure  of  merit,  critical  current  density  and  resistance,  when  the 
barriers  are  close  to  the  MIT;  and  (iv)  more  complicated  systems  like  SNSNS  junctions, 
and  so  on.  A  comprehensive  review  article  on  this  work  was  also  completed. 
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We  next  were  interested  in  determining  whether  there  was  any  universality  in  our  results 
that  could  be  practically  used  to  help  predict  properties  of  junctions  from  simple 
measurements  (perhaps  even  in  the  normal  state).  Motivated  by  quasiclassical 
calculations,  where  the  Thouless  energy  plays  a  prominent  role  (the  Thouless  energy  is  a 
quantum  mechanical  energy  derived  from  the  average  or  typical  time  that  an  electron 
spends  inside  the  barrier  of  the  JJ  and  Planck’s  constant).  In  the  quasiclassical  case,  there 
is  a  universal  formula  for  the  figure  of  merit  as  a  function  of  the  Thouless  energy.  When 
the  Thouless  energy  is  larger  than  the  superconducting  gap,  the  figure  of  merit  is 
proportional  to  the  gap,  while  in  cases  where  the  Thouless  energy  is  smaller  than  the  gap, 
it  determines  the  figure  of  merit.  We  found  that  all  of  our  metallic  barriers,  no  matter 
how  strong  the  scattering,  fit  the  quasiclassical  curve  perfectly.  When  we  went  through 
the  MIT,  we  found  a  simple  way  to  generalize  the  Thouless  energy  definition,  and  found 
the  generalized  result  also  fit  closely  to  the  quasiclassical  curve,  although  the  actual 
figure  of  merit  was  somewhat  higher  in  the  region  where  the  Thouless  energy  determined 
the  figure  of  merit.  We  feel  this  work  can  be  employed  both  for  diagnostic  and  design 
purposes  of  junctions.  To  get  the  best  figure  of  merit,  one  wants  to  keep  the  Thouless 
energy  larger  than  the  superconducting  gap,  if  at  all  possible. 

B.  Nonequilibrium  physics 

We  also  spent  a  significant  effort  on  developing  a  many-body  formalism  to  describe 
nonequilibrium  physics.  This  allows  for  an  understanding  of  the  nonlinear  response  of  a 
circuit  element  to  large  electric  fields  that  could  result  from  electronic  warfare,  or  natural 
sources  like  lightning.  Also,  since  nonlinearities  often  determine  the  limit  of  how  “hard” 
a  device  can  be  pushed  within  its  design  specs,  this  theory  also  allows  for  an 
understanding  of  how  devices  fail  under  ordinary  operating  conditions.  The 
nonequilibrium  formalism  was  developed  in  the  late  1950s  and  early  1960s,  but  to  date 
there  has  not  been  any  exact  solution  of  the  formalism  for  any  strongly  interacting 
system.  Instead,  much  work  has  focused  on  perturbative  effects,  which  may  not  be  valid 
if  the  interactions  are  strong  enough.  We  found  a  way  to  modify  equilibrium  dynamical 
mean  field  theory  to  work  in  the  nonequilibrium  case  as  well,  and  have  studied  two  of  the 
simplest  problems  with  these  techniques.  The  first  problem  is  a  “toy”  problem  that 
involves  determining  an  equilibrium  Green’s  function  from  the  nonequilibrium  approach, 
because  there  is  no  straightforward  way  to  directly  perform  the  analytic  continuation  to 
the  real  frequency  axis.  By  carefully  controlling  the  calculations,  via  exact  sum  rules,  we 
were  able  to  quantify  how  accurate  the  calculations  were,  and  areas  where  the 
computational  procedure  can  run  into  difficulties.  Armed  with  these  results,  we  next 
tackled  the  simplest  nonequilibrium  problem — the  nonlinear  response  of  a  strongly 
correlated  material  to  a  large  electric  field.  In  a  periodic  system,  the  noninteracting 
response  involves  the  so-called  Bloch  oscillations  (similar  in  many  respects  to  the  finite 
voltage  response  of  a  JJ,  where  a  dc  voltage  creates  an  ac  current).  As  we  turn 
interactions  on,  we  can  study  how  the  Bloch  oscillations  decay.  Although  a  simple 
Boltzmann  equation  approach  indicates  that  the  Bloch  oscillations  disappear  immediately 
after  scattering  is  turned  on,  our  data  indicates  that  the  oscillations  can  be  long  lived,  and 
maybe  even  survive  in  the  steady  state  if  the  field  is  large  enough. 


C.  Other  miscellaneous  work 


During  the  grant  period  two  comprehensive  review  articles  were  finished.  One  which 
summarized  how  to  use  inhomogeneous  dynamical  mean  field  theory  to  optimize  the 
performance  of  a  JJ  with  barriers  tuned  to  lie  close  to  the  MIT,  and  the  other  which 
summarized  over  a  decade  worth  of  work  on  DMFT  with  a  focus  on  the  Falicov-Kimball 
model.  In  addition,  a  book  is  being  written  on  the  theory  of  multilayered  nanostructures, 
including  JJs  and  should  be  complete  in  March  of  2006,  with  publication  in  the  winter  of 
2006-07. 

Other  work  included  development  of  theories  for  stripe  formation  and  phase  separation 
due  to  strong  electron-electron  interactions.  We  were  able  to  show  under  what  conditions 
phase  separation  or  stripe  formation  occurs.  This  is  of  interest  in  the  high-temperature 
superconductor  materials  which  may  have  stripe  formation  playing  a  role  in  the  materials 
properties.  We  investigated  spintronics,  by  developing  a  direct  approach  to  disordered 
problems  that  includes  the  spin-orbit  interaction  and  can  describe  the  Rashba  coupling. 

We  also  spent  a  significant  amount  of  time  on  thermal  transport  properties  in  the  hopes  of 
finding  insight  into  how  to  develop  low  temperature  thermoelectrics.  We  found  some 
promising  results  from  our  theory  when  one  is  close  to  particle-hole  symmetry,  but  has 
strong  correlations,  where  the  thermoelectric  figure  of  merit  can  get  large.  But  this  region 
requires  fine-tuning  of  parameters,  and  may  be  difficult  to  attain  in  real  systems.  We 
found  more  promising  results  at  high  temperature,  where  generically  power  generation 
appeared  to  be  viable  in  many  strongly  correlated  systems. 

Relevance  to  the  Superconducting  Electronics  program 

The  goal  of  this  work  was  to  provide  a  theoretical  framework  to  help  guide  experimental 
design  principles  for  optimized  JJ  performance.  Finding  a  new  technology  for  JJ  digital 
electronics  that  will  be  adopted  by  foundries  is  not  an  easy  task,  and  it  requires  a  dramatic 
improvement  of  performance  with  limited  extra  complication  in  the  processing  of  chips. 
Our  work  has  not  found  the  resolution  for  this  problem  to  yield  the  newest  JJ  for  digital 
electronics  chips.  Nevertheless,  our  work  has  shown  a  number  of  important  principles, 
which  include  (i)  when  the  barrier  is  on  the  metallic  side  of  the  MIT,  then  the  temperature 
dependence  of  the  figure  of  merit  is  the  most  important  quantity  for  determining  whether 
circuits  can  be  made.  Variations  in  the  critical  current  as  a  function  of  temperature  were 
quite  detrimental  for  the  TaxN-based  junctions  and  our  theory  indicated  that  this  is  an 
intrinsic  problem,  that  will  not  be  easily  overcome,  but  perhaps  could  be  addressed  by 
making  the  barriers  more  insulating,  (ii)  We  devised  a  new  metric,  based  on  a  generalized 
Thouless  energy,  to  use  in  determining  how  well  a  device  will  function.  Since  the 
Thouless  energy  can  be  determined  in  the  normal  state,  it  provides  a  simpler  diagnostic 
than  other  metrics,  but  it  has  not  yet  been  adopted  by  experimentalists.  Finally,  we 
worked  on  developing  a  new  formalism  that  can  shed  light  on  nonlinear  effects  in  many 
different  kinds  of  devices,  and  will  likely  find  wide  use  as  more  devices  with  strongly 
correlated  materials  are  made. 
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